value. A larger value of r’ will reduce 
the power dissipated in the load while 
increasing the reflection coefficient at 
the input port. It is now possible to ob- 
tain the excitation of the radiating 
waveguide closest to the load and the 
coefficients of the wave incident and re- 
flected at the input port of this cou- 
pling slot. The next coupling slot pa- 
rameter, r’, is chosen to realize the 
excitation of that radiating waveguide. 
One continues this process moving to- 
wards the source, until all the coupling 
slot parameters r’ and hence the Sn 
parameter of the 4-port coupler, r, are 
known for each coupling slot. The goal 


is to produce the desired array aperture 
distribution in the feed direction. From 
an interpolation of the computed mo- 
ment method data for the slot parame- 
ters, all the coupling slot tilt angles and 
lengths are obtained. From the excita- 
tions of the radiating waveguides com- 
puted from the coupling values, radiat- 
ing slot parameters may be obtained so 
as to attain the desired total normal- 
ized slot admittances. This process 
yields the radiating slot parameters, 
offsets, and lengths. The design is re- 
peated by choosing different values of 
r’ for the last coupling slot until the 
percentage of power dissipated in the 


load and the input reflection coeffi- 
cient values are satisfactory. 

Numerical results computed for the 
radiation pattern, the tilt angles and 
lengths of coupling slots, and excitation 
phases of the radiating waveguides, are 
presented for an array with uniform am- 
plitude excitation. The design process 
has been validated using computer sim- 
ulations. This design procedure is valid 
for non-uniform amplitude excitations 
as well. 

This work was done by Sembiam Ren- 
garajan of Caltech for NASA’s fet Propul- 
sion Laboratory. For more information, con- 
tact iaoffice@jpl.nasa.gov. NPO-48221 


© Time Manager Software for a Flight Processor 

John F. Kennedy Space Center, Florida 


Data analysis is a process of inspect- 
ing, cleaning, transforming, and mod- 
eling data to highlight useful informa- 
tion and suggest conclusions. Accurate 
timestamps and a timeline of vehicle 
events are needed to analyze flight 
data. When data is gathered onboard a 
rocket, precise time stamping is even 
more important due to the rocket’s 
high speeds and the requirement to in- 
tegrate data over time for inertial navi- 
gation calculations. 

Accurately time-tagging data cur- 
rently requires an additional accurate 
timecode generator board. This solu- 
tion is costly but is usually adequate for 
ground-based systems. However, this so- 
lution is not adequate for flight proces- 
sors on rockets. Rocket systems require 
more costly ruggedized equipment 
where weight, size, and power con- 
straints are an issue. Redundancy is also 
required, adding even more to the sys- 


tem’s weight, size, and power consump- 
tion. 

By moving the timekeeping to the 
flight processor, there is no longer a 
need for a redundant time source. If 
each flight processor is initially syn- 
chronized to GPS, they can freewheel 
and maintain a fairly accurate time 
throughout the flight with no addi- 
tional GPS time messages received. 
However, additional GPS time messages 
will ensure an even greater accuracy. 

Some modern microprocessors 
maintain a 64-bit internal time-base 
register that is incremented by a crystal 
oscillator, usually with a 20- to 100-MHz 
frequency. This time-base register can 
be read in an interrupt service routine 
(ISR) generated by the 1 pps signal 
from the GPS receiver. Next, a GPS 
time message is received. The time- 
base count is associated with the GPS 
time message time. 


When a timestamp is required, a get- 
time function is called that immedi- 
ately reads the time-base register. A 
delta count is calculated from the last 
GPS sync. The current time is calcu- 
lated by adding this delta time to the 
last sync time. This process calculates a 
timestamp with an accuracy measured 
in microseconds, depending on the 
processor clock speed and the accuracy 
of the processor clock. If a 1 pps GPS 
ISR is not available, the time base regis- 
ter can be synchronized with the re- 
ceipt of the GPS time message. If the 
microprocessor does not have a 64-bit 
internal time-base register, a count- 
down timer can be used. 

This work was done by Roger Zoerner of 
ASRC Aerospace Corporation for Kennedy 
Space Center. For more information, contact 
the Kennedy Space Center Innovative Part- 
nerships Office at 321-867-5033. KSC- 
13406 


© Simulation of Oxygen Disintegration and Mixing With 
Hydrogen or Helium at Supercritical Pressure 

These models can be used in simulations of gas turbine engines. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


The simulation of high-pressure tur- 
bulent flows, where the pressure, p, is 
larger than the critical value, p c , for the 
species under consideration, is relevant 
to a wide array of propulsion systems, 
e.g. gas turbine, diesel, and liquid rocket 


engines. Most turbulence models, how- 
ever, have been developed for atmos- 
pheric-/) turbulent flows. The difference 
between atmospheric-/) and supercriti- 
cal-/) turbulence is that, in the former sit- 
uation, the coupling between dynamics 


and thermodynamics is moderate to 
negligible, but for the latter it is very sig- 
nificant, and can dominate the flow 
characteristics. The reason for this stems 
from the mathematical form of the 
equation of state (EOS), which is the 
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perfect-gas EOS in the former case, and 
the real-gas EOS in the latter case. 

For flows at supercritical pressure, p, 
the large eddy simulation (LES) equa- 
tions consist of the differential conserva- 
tion equations coupled with a real-gas 
EOS. The equations use transport prop- 
erties that depend on the thermody- 
namic variables. Compared to previous 
LES models, the differential equations 
contain not only the subgrid scale (SGS) 
fluxes, but also new SGS terms, each de- 
noted as a “correction.” These addi- 
tional terms, typically assumed null for 
atmospheric pressure flows, stem from 
filtering the differential governing equa- 
tions, and represent differences between 
a filtered term and the same term com- 
puted as a function of the filtered flow 
field. In particular, the energy equation 


contains a heat-flux correction (q-cor- 
rection) that is the difference between 
the filtered divergence of the heat flux 
and the divergence of the heat flux com- 
puted as a function of the filtered flow 
field. In a previous study, there was only 
partial success in modeling the q-correc- 
tion term, but in this innovation, success 
has been achieved by using a different 
modeling approach. 

This analysis, based on a temporal 
mixing layer Direct Numerical Simula- 
tion database, shows that the focus in 
modeling the q-correction should be on 
reconstructing the primitive variable 
gradients rather than their coefficients, 
and proposes the approximate deconvo- 
lution model (ADM) as an effective 
means of flow field reconstruction for 
LES heat flux calculation. Further, re- 


sults for a study conducted for temporal 
mixing layers initially containing oxygen 
in the lower stream, and hydrogen or 
helium in the upper stream, show that, 
for any LES, including SGS-flux models 
(constant-coefficient Gradient or Scale- 
Similarity models, dynamic-coefficient 
Smagorinsky/Yoshizawa or mixed 
Smagorinsky/Yoshizawa/ Gradient mod- 
els) , the inclusion of the q-correction in 
the LES leads to the theoretical maxi- 
mum reduction of the SGS heat-flux dif- 
ference. The remaining error in model- 
ing this new subgrid term is thus 
irreducible. 

This work was done by Josette Bellan and 
Ezgi Taskinoglu of Caltech for NASA’s fet 
Propulsion Laboratory. For more informa- 
tion, contact iaoffice@jpl.nasa.gov. NPO- 
47645 
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